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Abstract 

The Inverse Amplitude Method is a powerful unitarization technique to enlarge the energy appli- 
cability region of Effective Lagrangians. It has been widely used to describe resonances from Chiral 
Perturbation Theory as well as for the Strongly Interacting Symmetry Breaking Sector. In this work 
we show how it can be slightly modified to account also for the sub-threshold region, incorporating cor- 
rectly the Adler zeros required by chiral symmetry and eliminating spurious poles. These improvements 
produce negligible effects on the physical region. 

1 Introduction 

Effective Field Theories provide a systematic and model independent approach to systems whose symmetries 
and low energy degrees of freedom are known but whose description in terms of an underlying fundamental 
Quantum Field Theory is out of reach. The two cases of interest for this work are, on the one hand, Chiral 
Perturbation Theory (ChPT) [1] that describes effectively the low energy dynamics of hadrons, unaccessible 
to perturbative QCD calculations in terms of quarks and gluons and, on the other hand, the effective 
description of the Strongly Interacting Electroweak Symmetry Breaking Sector (SISBS) [2], whose underlying 
fundamental theory remains unknown. 

Both cases have in common the existence of an spontaneous symmetry breaking of a global chiral 
SU(N)l x SU(N)r group down to an SU(N)l + r group. The Goldstone Theorem implies the presence 
of N 2 — 1 massless Goldstone Bosons (GB) in the particle spectrum. These GB thus become the relevant 
degrees of freedom of the system below a chiral scale A x , where an Effective Chiral Lagrangian can be built 
in terms of just those GB as the most general derivative expansion consistent with the known symmetries. 
Note that there is a well defined power counting so that divergences generated in loop diagrams that use 
vertices up to a given order can be renormalized into the coefficients of higher order terms. In this sense 
these Effective Theories are renormalizable order by order. 

In addition, in both examples, there are small explicit symmetry breaking terms. For QCD the relatively 
small masses of the lightest three quarks provide a mass M « A x to the GB, identified with the pions, 
kaons and etas, so that the effective approach becomes in practice a derivative and mass expansion. For 
the SISBS, there is a local SU(2)l x U(l) symmetry whose gauge bosons couple to the '"would-be GB"' 
that, in a unitary gauge, disappear from the spectrum, giving rise to gauge boson longitudinal components 
that thus acquire a mass My. In this way the SU(2)l x U(l) gauge symmetry is spontaneously, but not 
explicitly, broken to the electromagnetic group U(l) em- In this case, in addition to the derivative expansion, 
one expands also in terms of electroweak coupling constants g and g' . The so-called Equivalence Theorem 
(ET) [3] states that at high energies (in gauges, and to leading order in momenta over My and g and 
<?') amplitudes involving longitudinal gauge bosons can be calculated as if they were GB, which, being 
pseudoscalars, are much easier to handle. Although this is a high energy limit, there is a generalization to 
the Effective Lagrangian formalism [4], that, for practical purposes, allows us to identify, up to the difference 
in scales, the formalisms of SU(2) ChPT and the SISBS and therefore, from now on we will be referring to 
ChPT, but keeping in mind that our results have a straightforward translation to the SISBS. 
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Both cases above are examples of strongly interacting systems whose most salient feature is the saturation 
of unitarity and the associated resonant states, which lie beyond the reach of perturbative energy expansions. 
Thus, it may seem that the use of effective Lagrangians is limited to energies below those resonances, whose 
effects are encoded in the values of higher order effective coupling constants. However, since unitarity fixes 
the imaginary part of inverse partial waves in the elastic region, the effective Lagrangian approach is also 
useful in the resonant region, for instance, used inside a dispersion relation, in order to obtain the rest of 
the amplitude. These techniques are known as unitarization methods, and reproduce simultaneously the low 
energy expansion and the lightest resonances without including them explicitly in the Lagrangian. The great 
advantage is that such resonances and their properties are generated without prejudices about their nature 
or their existence. Also, since the Lagrangian symmetries and some features of the effective constants can 
be directly related to the underlying theory, like QCD, one can study the properties of these states based 
on more fundamental grounds. One of the most extensively used unitarization techniques is the Inverse 
Amplitude Method (IAM) [6] [7l [8] , that uses the fully renormalized effective chiral expansion, without 
any further approximation and without introducing any other spurious parameter, but just the effective 
constants up to a given order. Within hadronic physics, it generates the well known vector resonances and 
the more controversial scalars using parameters consistent with one loop ChPT, allowing to establish their 
different nature in terms of their dependence on the number of colors. It has also been possible to extend it 
up to two loops [9] , the finite temperature formalism [TUl [11] and to the pion-nucleon sector [12] . Within 
the SISBS [13 , it provides the prediction of the general resonance spectrum and how well it could be detected 
at LHC. 

However, it is known pj Q3] that the IAM fails to reproduce correctly the Adler zeros that appear in 
the subthreshold region of some partial waves as a consequence of chiral symmetry [TS]. Furthermore, it 
generates spurious poles, or "ghosts" , thus questioning its reliability in that region, and also casting some 
doubts about the robustness of the results in the physical region if such structures were properly accounted 
for. 

The aim of this paper is to show that a very simple modification of the IAM can take into account correctly 
those zeros and ghosts. This modification corresponds to terms that had been neglected in the original 
dispersive derivation of the IAM since they contribute to higher orders in the chiral expansion. Actually, we 
will check explicitly that such procedure is justified in the physical region, where these modifications yield 
negligible contributions, thus showing the robustness of the standard IAM. However, apart from improving 
the IAM consistency, these terms are essential in the sub-threshold region which is relevant to study the 
effect of chiral symmetry restoration on resonances |11) , or their dependence on quark masses |16j . 

In the next section we will thus revisit the standard IAM derivation from dispersive theory, where Adler 
zeros are neglected, paying special attention to the role of those zeros in the subthreshold region. In section 
[3]we will present a naive way of extending the IAM amplitude, without using dispersion relations, that solves 
the caveats in that region. Section [4] will show a dispersive derivation of a more general modified amplitude, 
for the case of equal masses (e.g. hit scattering). The case of unequal masses, like in ttK scattering, deserves 
a separate discussion, for the reasons explained in section [5] Finally, in section [5] we will present some 
numerical results for the modified amplitudes. 

2 The Inverse Amplitude Method 

2.1 Dispersive derivation 

The one-channel Inverse Amplitude Method (IAM) [5][5][71[§] can be obtained by using ChPT up to a given 
order inside a dispersion relation. To simplify the discussion, let us first consider pion-pion scattering partial 
wave amplitudes of definite isospin / and angular momentum J, although for brevity we will simply call 
them t, whose analytic structure in the s-plane is shown in Fig. [1] The physical right hand cut comes from 
unitarity and starts at threshold Sth, while the left hand cut comes from the t,u channels. 

The inverse of t(s) has the same analytic structure, except for the possible presence of poles corresponding 
to zeros of t(s). In particular, chiral symmetry requires the existence of the so-called Adler zeros below 
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Figure 1: Pion-pion scattering partial waves analytic structure and integration contour. 



threshold in scalar partial waves [UHI], the case of interest in this work, which we denote by sa- Hence it 
is then possible to write the following dispersion relation for the inverse amplitude: 

1 s-zq r d Im lA(z) +Lcm + pc{m (1) 



t( S ) t(z ) TT J Sth (Z-S)(Z-Z ) 

here and in the following, we will simply write z instead of z + ie with e > for the imaginary parts inside 
the cut integrals. Note that we have explicitly written the integral over the right hand cut (or physical cut, 
extending from threshold, s t h to infinity) but we have shortened by LC the equivalent expression for the left 
cut (from — oo to 0) and the pole contribution. We could do similarly with other cuts, if present, as in the 
ttK case. In addition we have made one subtraction to ensure convergence, at a point Zq ^ sa, s 2 - 
We now recall that unitarity, for physical values of s in the elastic region implies: 

Imi(s) = a\t(s)\ 2 => lm-^ = -a(s), (2) 

where er(s) = 2pcM / \fs. Let us remark that since Iml/i = —a we know exactly the integrand over the 
elastic cut. 

In contrast, ChPT amplitudes are obtained as a series expansion t(s) = £2(5) + t^s) + ... where t 2 (s) — 
0(p 2 ), t^s) — 0(p 4 ) and p stands for the pion mass or momentum. Therefore elastic unitarity is not satisfied 
exactly, but only order by order as follows: 

lmt 2 (s)=0, Imt 4 (s) =a(s)\t 2 (s)\ 2 , ... (3) 

Let us also note that £2(5) is a pure polynomial and has no cuts, and we can thus write a trivial dispersion 
relation for l/£ 2 (s) that reads 

where now the pole contribution is due to S2, the Adler zero of £2- In addition, except for the poles, the 
function t4(s)/t?,(s) has the same analytic cut structure of 1/t and, using Eqs.© and Q, over the physical 
cut we find 

*4(s) i , T 1 
9 . . = <J(s) = — 1m — - . 
t 2 2 (s) t(s) 

We can therefore write another dispersion relation similar to that of l/t(s), but for t4(s)/t?,(s) 



l m -^-L =(7 ( s ) = -l m _. (5) 



t 4 (s) U(z ) , s-z r & Enfr(*)/ta(*)» +LC{u/tl)+PC{u/t l l (6) 



t 2 (s) 2 t 2 {z ) 2 7T J s (Z-S)(Z-Z ) 
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where the pole contribution, once again, is due to the Adler zero of t%. 

We are now going to relate the dispersion relation for l/t(s) with that for ti(s) /t\{s). As we already 
commented Iml/t(s) = — Imt^/t^s) on the right cut, and therefore the integrals over the physical cuts for 
l/i(s) and ti/t\(s) are exactly opposite to each other. In addition, using ChPT we find that LC{\/t) ~ 
-LC{U/tl), which is a well justified approximation, since, due to the subtraction, the integrand of LC 
is weighted at low energies, precisely where ChPT applies. Finally, we have to evaluate the subtraction 
constant in Eq.(fT]), and this can only be done as long as zq is in the low energy region, where it is perfectly 
justified to use ChPT to find l/t(zo) ~ l/t2(zo) — t4(zo)/t2{zo) 2 . However, note that this expansion is a very 
bad approximation for zq near S2 or sa, where t 2 and t vanish. Therefore, we only know how to relate those 
dispersion relations for subtraction points zq in the low energy region, but sufficiently far from the Adler 
zeros. In section 6 we will check that the results have very little sensitivity to the choice of zo as long as it 
lies in this region. When this is the case, using Eqs.JI]), (TJ| and ([6]) we can write 1/t as 

-L ~ J- _ p& _ pc(i/t 2 ) + pc(u/4) + pc(i/t) (?) 
t(s) t 2 {s) t 2 {sy 

The standard dispersive derivation of the IAM [6l [7] simply neglected the sum of pole contributions to 
arrive at 

t IAM (s) ^ (8) 

thus providing an elastic amplitude that satisfies unitarity and has the correct low energy expansion of ChPT 
up to the order we have used. When such amplitude is chirally expanded to 0(p 4 ), this implies [7] that the 
total pole contribution to t(s), even without its explicit calculation, has to be 0(p 6 ). In sections I3I4I5I we 
will calculate it explicitly to arrive to the modified IAM. 

2.2 IAM properties and its naive derivation 

Incidentally, we can recast Eq.J^J) as: 

t(s) = 5 -17 \ —FT , (9) 
Ret i (s) — lo-(s) 

and thus it seems that the IAM can also be derived in a much simpler way by replacing Rei -1 by its 0(p 4 ) 
ChPT expansion Ret' 1 = (t 2 - Rei 4 )/i|- This is the way unitarization methods are usually presented, 
although it makes no use of the strong analytic constraints of amplitudes, which are indeed absent in Eq.([9j). 
Furthermore, immediately the criticism is raised that the ChPT expansion cannot be used at high energies. 

However, note that the dispersive one-channel IAM derivation in the previous section imposes analyticity 
in the form of two dispersive integrals and makes use of the ChPT expansion (up to one-loop in this case) for 
the subtraction constants and the left cut. The use of ChPT for the subtraction constants is well justified, 
since it is used at s = zq in the low-energy region. Since the integral extends to infinity, ChPT may seem a 
worse approximation for the left cut or possible inelastic cuts, however, the subtraction suppresses the high 
energy region that contributes little, as explained above. Furthermore, when the IAM is used for physical 
values of s above the physical threshold, the left cut is damped by an additional l/(z — s) factor, and not only 
the high energy part, but its whole contribution is rather small. There are no model dependent assumptions, 
but just approximations to a given order, and therefore the approach provides an elastic amplitude that 
satisfies unitarity and has the correct ChPT expansion up to that given order. It also is straightforward to 
extend it to higher orders [7J [5] ■ 

Moreover Eq.([2]) is only valid on the physical cut, whereas the dispersive derivation allows us to consider 
the amplitude in the complex plane, and, in particular, look for poles of the associated resonances. Actually, 
already ten years ago [7] the poles for the p(770), if* (892) and most interestingly, the controversial a (also 
called /o(600)), where generated without any model dependent assumptions. 

Obviously, and contrary to a wide belief, the IAM contains a left cut and respects crossing symmetry up 
to, of course, the order in the ChPT expansion that has been used. The confusion may come from the fact 
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that the IAM has also been applied in a coupled channel formalism, for which there is still no dispersive 
derivation and has been frequently used by approximating further the amplitudes neglecting tadpole and 
left cut terms |T7]. But, strictly speaking that would not be the full IAM, which definitely has a left cut. 

2.3 IAM caveats: Adler zeros and ghosts 

To end this section, we recall that in the dispersive derivation the sum of pole terms (PC) in Eq.([7} is 
neglected, since it yields a higher order contribution [7]. However, this leads to a couple of problems related 
to the presence of the Adler zeros below threshold in the scalar waves. Let us first note that despite we have 
used in the IAM the ChPT expansion up to next to leading order (NLO), due to the t 2 (s) 2 factor in the 
numerator, it vanishes at s 2 , which is only the leading order (LO) chiral approximation to the exact Adler 
zero sa- In addition it is a double zero instead of a single zero. This, as we will see, leads to the appearance 
of a spurious pole on the real axis close to the Adler zero. As a consequence, the predictions of the standard 
IAM below threshold and in particular, around the Adler zero position, are not reliable. 

In fact, note that since the interval s £ (0, Sth) lies within the convergence region of the chiral expansion 
and t 2 changes sign at s = s 2l t 2 — t^ turns out to be, for the cases of interest here, a continuous, monotonically 
increasing (or decreasing) function in (0, Sth) that changes sign from s = to s — sth- Therefore, there is 
one single point s where the denominator of j8]) vanishes, which, as long as £4(^2) 7^ 0, produces the spurious 
pole below threshold discussed above, i.e., a not-observed tttt bound state. 

For instance, let us consider I = J = ttit scattering (for the values of the low-energy constants quoted 
in section [6l This is an attractive channel, so that ti — t± is positive at threshold. Since that function 
is negative at s = and at s — s 2 , s > s 2 in that case, as showed in the upper right panel of FigJ^l we 
find s ~ (110 MeV) 2 and s 2 — (99 MeV) 2 . In the / = 2, J = channel, which is repulsive, t 2 — £4 is 
positive at s = and negative at s = s 2 (^4(^2) > 0) so that s < s 2 . In that channel, s ~ (196 MeV) 2 and 
s 2 ~ (198 MeV) 2 as seen also in FigH] 

Let us also point out that, together with the first Riemann sheet spurious pole just discussed, the IAM 
has a companion pole in the second Riemann sheet below threshold. For instance, the second-sheet IAM for 
7T7T scattering reads: 

±( \IAM,IIsheet ^(^) 2 

W t 2 {s)-U{s)-2d{s)tl{s) ( ' 

where <r(s) = icr(s — i0 + ) = ■\/4m 2 r /s — 1 for < s < 4m 2 . Thus, if we are dealing with an attractive 
channel, like the 00 one, the denominator of (fl"0|) is positive near threshold (dominated by t 2 > 0) and 
diverges to minus infinity as s — ► + , so that it must have at least one zero, which again generates a pole if 
£4(52) 7^ 0. Since — at 2 in Eq. (fTU|) is also an increasing function, there will be only one such pole. 

In conclusion, the existence of the perturbative Adler zeros and the fact that the IAM amplitude does 
not reproduce them well leads to the presence of spurious poles. Similar conclusion had been noticed in 
[18j . In the next sections, we present, first, a very simple construction of a modified IAM, along the lines of 
the previously discussed naive derivation of the standard IAM obtained without using dispersion relations, 
that solves these problems and next we show two dispersive derivations of the modified IAM. One of them 
is subtracted at arbitrary zq (within the range of validity of our approximations) and the other one at 
the Adler zero. We will show how the modified IAM obtained naively corresponds to a particular limit 
of the first dispersive approach and comes out directly in the second. Finally, we will also show that the 
differences between the modified IAM formulae are negligible numerically, and that, while fixing the above 
mentioned problems, the modified IAM does not yield any significant modification over the standard IAM in 
the physical region and to the resonance poles. Therefore, the results obtained so far in the literature with 
the IAM remain valid. 
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3 Modified IAM: Naive derivation 



First of all we will set some notation: as before, we denote by sa the Adler zero of the "complete" partial 
wave, namely t(sA) = 0. In addition, we will also use the approximations to the Adler zero at LO, s 2 , and 
NLO, s 2 + S4- Thus, £2(^2) = and t 2 (s 2 + s 4 ) + t 4 (s 2 + s 4 ) = 0. 

In this section we present a naive, and intuitive, derivation leading to a partial wave definition which 
does not have the Adler zero related problems discussed above. The derivation follows closely [TJJ, where 
this problem was addressed in the context of real pion scattering poles arising by medium effects such as 
temperature and density. In that paper, there was not a formal proof based on dispersion relations, as the 
one we will present here later, and it was limited to pion-pion scattering. 

/From the discussion in the previous sections, it is clear that if we modify the inverse amplitude as 
\/t IAM (s) — > l/t IAM (s) + A(s)/t 2 with A(s) an analytic function at least off the real axis, real for real s, 
the unitarity and analytic properties of the amplitude remain unaltered. The modified IAM, (from now on 
called mlAM), reads then: 

t mIAM( \ = i(f) (n \ 

{) t 2 (s)-U(s)+A™i*M( s y [ 1 

Consider now the case of mr scattering, where we have simply t 2 (s) — t 2 (s 2 )(s — s 2 ) with t' 2 {s 2 ) constant. 
Then the Laurent expansion around s = s 2 of the standard IAM reads: 



t IAM (s) - ^( S2 ) 2 ( S _ S2 )2 t> 2 { S2 f{ S -S 2 ) 

The idea is that if we want the amplitude to have only an Adler zero of order one at s = sa, we must 
subtract from l/t IAM the above double and single pole contributions at s — s 2 and add a single pole at sa, 
i.e, 

1 1 U{s 2 ) t' 2 (s 2 ) - t' 4 (s 2 ) c 



t mIAM {s) t IAM {s) ^ (sa)2(s _ S2)2 t / (s2 )2 (s „ S2) s _ SA < 



(13) 



where c is a so far undetermined constant that, as we will show now, can be fixed by demanding that the 
mlAM formula matches the perturbative ChPT series to fourth order, namely A = 0(p 6 ). 

In practice, it is simpler to keep track of the different chiral powers by counting the powers of /~ 2 , where 
/ is the pion decay constant. Thus, since t' 2 (s 2 ) = 0(f~ 2 ), t 4 = 0(/ -4 ), and s 4 — 0(/ -2 ), expanding the 
expression t 2 (s 2 + s 4 ) + t 4 (s 2 + s 4 ) — around s 2 we find: 

s 4 = -t 4 (s 2 )/t' 2 ( S2 )+0(f- 4 )- (14) 

Using this in Eq. (fT3"]) with 5^4 = 52 + 54 + 0(/ -4 ), and requiring that Eq. fTTj) matches the chiral expansion 
at low energies, we find the first two orders of the chiral expansion for c: 

which leads to 

AmIAM/ \ . I \ ( s 2 — Sa){s — S 2 ) , . . . 

A {S)=t 4 (s 2 ) [t 2 ( s 2) - t 4 (S 2 ) . (iDj 

S - SA 

Therefore, the mlAM in Eq. fTTj) with A mIAM (s) in Eq. (fT6]) and sa approximated by its chiral expansion 
given above, matches the chiral expansion of the amplitude up to fourth order and has the Adler zero at 
the same position and with the same order as the perturbative amplitude. Furthermore, we have solved in 
turn the spurious pole problem. In fact, since A mIAM (s 2 ) — t 4 (s 2 ), the denominator of Eq. fTlj) vanishes at 
s = s 2 . But, for s s 2 we have shown that A mI (s) = 0(f~ 6 ) and therefore our previous argument about 
the monotonous behavior of the denominator still holds, so that the denominator vanishes only at s = s 2l 
where there is no pole contribution. The same holds for the spurious second-sheet poles. In Fig. [2] we show 
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the mlAM amplitude in the I = J = channel, which we observe that is not singular below threshold and 
remains close to the standard IAM result away from the Adler zero region. The same situation takes place 
in the / = 2, J — channel. Finally, we have checked that, as expected from our previous arguments, the 
Jo (600) or a pole remains at the same place either using the second Riemann sheet extensions of the mlAM 
or the IAM (see Tabled]). 

In the next sections, we will check how a modified IAM can also be obtained by considering explicitly the 
pole contributions in the dispersive derivation, thus ensuring the correct analytic properties of the amplitude. 
We will also show that the modified formula obtained with the naive derivation in this section can also be 
obtained as a particular case. 



4 Modified IAM: dispersive derivation for equal masses. 
4.1 Pole Contribution to the standard derivation 

The derivation of the modified IAM from dispersion theory follows that in section 2 up to Eq.([7]), but keeping 
the pole contributions, which, by evaluating the corresponding residues, read: 

PC(l/t2) = TH-T ( T 1 - - (17) 



t' 2 {s 2 ) \s - s 2 z Q - s 2 

pew® = rH(r^-r^) + rH(- -) < 18 » 

PC(l/t) = -J— A— —) (19) 



t'(s A ) \s- s A zo- SA 

where we have assumed a single zero in t(s^), as the presence of the l/(s — Sa) factor in PC(l/t) shows. 
As we have discussed above, and as it is detailed in the Appendix, by expanding chirally —PC(l/t 2 ) + 
PC(t4/t 2 ) + PC(l/t), the poles contribute to t(s) at higher order, and that is why they were customarily 
neglected. However, these pole contributions contain the terms needed to have the Adler zero in the correct 
position. In addition, by taking the limit s — > s 2l (~PC(l/t 2 ) + PC^/if)) tends to the term needed to 
cancel the double zero of the IAM in s 2 , see Eq. f(T2|) . and the spurious pole will also disappear. 

In summary, the modified IAM obtained from dispersive relations subtracted at zq, that we will denote 
by zqIAM, can be written again as: 

.z IAM, n ^|(f) (2Q) 



where now 

A ZoIAM (s) = A mIAM (s) ~ h^f A mIAM (z ) (21 ) 



Of course, the position of the Adler zero for the total amplitude is not known, but since we have been working 
with ChPT to one loop, we can impose the Adler zero to be located in its one-loop position, namely, in the 
above formulae, we have to replace sa — * s 2 +s 4 which is obtained from the equation t 2 (s 2 +S4)+t4(s 2 +s 4: ) = 
as explained in Section [31 Thus, to obtain (|2"Tj) we have made use of t 2 (s) — t' 2 (s 2 )(s — s 2 ) and chirally 
expanded l/t'(sA) — l/t' 2 (s 2 ) — t' 4 (s 2 ) / t' 2 (s 2 ) 2 , which is perfectly justified near sa- Note that we can use the 
chiral expansions around sa and s 2 because we no longer expand the inverse of the amplitudes but that of 
their derivatives, as they appear in the residues of the pole contributions. 

Note that, once again, the factor A mIAM (s) that appeared in the previous section is present, but now 
there is an additional and very similar piece that carries a zq dependence, which occurs due to our truncation 
of the ChPT series when approximating the subtraction constants and pole contributions. This additional 
term in A Z ° IAM is also 0(/~ 6 ) but, as we will see below, as long as zq lies within the range where our 
approximations remain valid, it is numerically small not only in the physical region but also below threshold. 
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Therefore, this term can be dropped without spoiling the right chiral behavior of the amplitude in the 
subthreshold region, obtaining again the mlAM in section [3l which, for the moment, we have justified only 
numerically. 

Furthermore, it is tempting to take the z$ — > oo limit in Eq. (12"T]) and recover the mlAM of the previous 
section by noting that the second term in Eq. (|2ip vanishes in that limit. However, this is just a formal 
justification of our naive derivation since, actually, we required for the subtraction point zq to remain in 
the low energy applicability region of ChPT. Nevertheless, we will see in the next section that there is an 
alternative dispersive derivation of the mlAM, which is somewhat different from the standard dispersive 
derivation since it requires subtractions at the Adler zero of each function. 



4.2 Subtraction at the Adler zeros 

The way to derive Eq. (Til)]) from dispersive theory is to make the subtractions at the place where we already 
have a pole. Note, however, that l/t has its pole in the Adler zero at s A whereas t\jt\ has its pole at s 2 , so 
that we have to write 

' - S ~ SA f dz 1^1 -+LC(l/t) + PC(l/t), (22) 



t(s) 7T J RC (z - S A )(z - s) 

U(s) _s-s 2 f o(z) + LC {t i /4) + PC(U/4), (23) 



t 2 (s) 2 7T J RC (z~S 2 )(z~s) 

where, for brevity, we have already used the clastic unitarity condition Eq. ([5)1 . As usual, we will approximate 
LC(l/t) ~ —LC(t4,/t 2 ), since that is the result of the chiral expansion at low energies where the integral is 
weighted. In the above relations, the pole contributions PC (l/t) and PC(ti/t%) come now from a double 
and triple pole respectively, and read 

PC (l/t) = ,, J r- *", (aA L (24) 

V 1 1 t'(s A )(s-s A ) 2t>( SA y , 1 ' 

PC(ti/t\) — — -2 + -j-. 4 J r + —77 rj- (25) 

t' 2 (s2r(s - s 2 y t' 2 (s 2 y(s - s 2 ) 2t' 2 (s 2 y 

In addition, it is important to remark that these pole contributions diverge either on s 2 or s A , so that, for 
s ~ s 2 or s ~ s A they are, by far, the dominant contributions, since at the same time the right and left cut 
terms tend to zero, and therefore become negligible. 

Outside that region, the other terms become relevant and the difference with our previous derivations 
is that now we also approximate the l/t integral term over the right cut by using (s — s A )/(z — s A ) ~ 
(s — s 2 )/(z — s 2 ), which is its LO chiral expansion. This is a remarkably good approximation for the 
dispersion relation as long as z is sufficiently far from s 2 and s A , which is indeed the case right cut integral. 
Of course, the l/t right cut term should vanish at s A and now it does not, but, as we have just commented, 
the pole contribution diverges precisely at s A and thus is largely dominant over the integral, which therefore 
can be completely neglected. 

Once again, we simply add Eqs. ([22|) and (|23| to obtain the mlAM equation: 

1 U(s) , ti(s 2 ) , t' 4 (s 2 ) , 1 



V nIAM {s) t2(s) 2 t ' 2 (s 2 f(s-S 2 f t' 2 (s 2 f( S - S 2 ) t'( SA )(s-S A Y 

where, in the pole contributions, we have used that 

-*M + «W =0(r»), (27) 

2t'(s A y 2t' 2 (s 2 y 

which, once again, can be safely neglected since it correspond to the chiral expansion at very low energies. 
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Finally, if we evaluate perturbatively 



F(b"^y"l^ +0(r2) ' (28) 

and add to Eq. (j2l))) written as l/t 2 (s) — l/(^ 2 (s 2 )(s — s 2)), we obtain for t(s) Eq. ([TT|) with A(s) given by 
Eq. tflTj)) . i.e, we recover our naive derivation of the modified IAM by subtracting at the Adler zeros. 

Note however, that in contrast to the derivation where we used an arbitrary zo, now all subtractions 
have been performed in the very low energy region so that chiral expansions for pole terms are well justified 
within ChPT. The price to pay is that in the right cut integral terms we have approximated sa by s 2 , which 
is irrelevant in the Adler zero region since the pole contributions dominate there, and a remarkably good 
approximation in the physical and resonance regions. Of course, we have still used exact elastic unitarity in 
the integrands over the physical cut, which ensures that the modified IAM satisfies exact elastic unitarity. 



5 Modified IAM: unequal masses. 

When dealing with unequal masses, as in the irK scattering case that we will use for reference, in addition to 
the left and right cuts, there is also a circular cut centered at s = with radius y m 2 K — m\ that contributes 
to the dispersion relation. This circular cut lies in the low energy region within the applicability range of 
ChPT and thus, as we do with the left cut, we will approximate the inverse amplitude by its ChPT series 
to fourth order. Taking into account that t 2 (s) has no cuts we will have CC(l/i) = —CC{t^/t%), where CC 
stands for the circular cut contribution. This is the same approximation used before for the left cut and 
therefore we obtain the same IAM dispersive derivation. Hence, there is still the same problem with Adler 
zeros and spurious poles. 

The solution given in previous sections works similarly well for the 7 = 3/2 and J = partial wave. 
However, for the I = 1/2, J = irK scattering channel complications arise due to the form of the LO partial 
wave t 2 {s) which has two zeros instead of one. 

t>(s) = - 5iS - S ^ iS - S2 -\ (29) 

128 7T/S 

with s 2 ± = | ( m K + m Tc ^ m K ~ 7m^-m| + 4m£j whose values are s 2+ = 0.24 GeV 2 , and s 2 _ = 

—0.13 GeV 2 . In particular, this means that, contrary to the previous cases t"(s 2 ±) 7^ 0, which complicates 
the derivation of a modified IAM. Nevertheless, once again we have found for this special case of the I = 1/2, 
J = channel, a naive and two dispersive derivations, that we detail next. 

5.1 Naive derivation, 1 = 1/2, J = channel 

Following section 3, we define \jt mIAM = l/t IAM + A mIAM (s)/4(s) with A mIAM (s) an analytic function 
at least off the real axis and real for real s to preserve the unitarity and analytic properties of the original 
amplitude. Next we expand l/t IAM in Laurent series around s = s 2+ , taking into account that now £ 2 (s) 7^ 0, 
obtaining 

1 _ hM , t' 2 (s 2 +) - t' 4 (s 2 +) , t 4 (s 2+ )£ 2 (s 2+ ) , O ( s _ S2+ )0 ( 30 ) 



t^(s) t' 2 { S2+ f{s s 2+ f t' 2 { S2+ Y{ S ~ S2+ ) t' 2 (s 2+ r( S ~s 2+ ) 

As in section 3 we subtract the pole at s 2+ and add a pole at sa to the inverse amplitude, and the constant 
in the sa pole term is calculated perturbatively in order to match the chiral expansion at low energies. Note 
that, following our previous arguments, we do not need to subtract the s 2 _ pole, since Imt^ 7^ on the LC, 
so that no spurious pole will appear in that region. Proceeding then as in section [3l we get now: 



1 ml AM 



(s) 



t' 2 ( S 2 + ) 2 



U(s 2+ ) _ 
(s-s 2+ ) 2 - s 2+ )(s - s A ) 



{s 2+ - s A ) ( ,. , . . U(s 2+ )t 2 (s 2+ ) \ 

r W2+ - t i( s 2+) + 777 \ 

-s 2+ ){s-s A ) V t' 2 (s 2 +) J 



(31) 
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Then, Eq. fTTj) with A mIAM (s) above unitarizes the I = 1/2, ,7 = channel and has an Adler zero of the 
correct chiral order and has no spurious pole. Let us remark that the above A mIAM (s) is a generalization 
also valid for the equal mass case, since it is reduced to the A mIAM (s) we already obtained in Eq. p^)) when 
%(s 2+ ) = 0. 

Apart from this naive formal derivation we can follow a dispersive approach, detailed next, in which we 
make use of analyticity and the ChPT series is only used within its applicability region. 



5.2 Dispersive derivation at z , I = 1/2, J = channel 

The derivation follows exactly that of section 2.1 but now, due to the additional zero, s 2 _, of t 2 , which 
lies on the negative axis, we cannot simply write, as usually done, Disc £4/^ = 2i j LC lmt4,/t 2 , nor 
J LC Disc l/<2 = 0, where Disc / = f(x + ie) — f(x — ie) with real x. This zero also modifies the form of the 
dispersion relation for l/t2, that now reads 

' ' + PC + (l/t 2 )+PCL(l/t 2 ), PC^i/ia) = —U- (_i 1 ) (32) 



t 2 (s) t 2 (z ) '' ' ■ " ' " t' 2 (s 2 ±) V S - S 2± Z - S 2 ± 

Following the same steps as in section 2.1, we approximate 1/t on the left cut by its chiral expansion, but 
now taking into account that, as mentioned above, J LC dzDisc (l/t 2 (z)) 0, i. e. 

LC(l/t) = S ~ Z ° f Disc V^Cg) ^ s ~ z o f Disc(lA 2 (z) -U{z)/t 2 {zf) 
2ni J LC (z - s)(z - z ) ~ 2ni J LC (z - s)(z - z ) 



where we now get —LC(t4/t 2 ) as before, and a new term coming from l/t 2 , 
s — z f Discl/t 2 (z) s — z f (l/*2 (z + ie) — l/t 2 (z — i e) ) 



27r« Jlc i z - s )( z - z o) 2?ri J LC (z - s)(z - z ) 



(34) 



f (z-s 2 -)/t 2 (z) 

= -{.S-Zo) r- 5(Z - S 2 -) - , , 

Jlc \ z s )\ z z o) * 2 ( s 2-) V s ~ s 2- z o - s 2 - 

which is equal to PC-(l/t 2 ), so they will cancel in the expression for 1/t. Thus we again obtain Eq.([7j), but 
now the explicit expression for PCfa/t^) at s = s 2+ is different from Eq. p^|) because t 2 (s) ^ 0, 

, 2 \ _ t 4 (s 2 +) ( 1 1 \ , f f 4( s 2+) t A (s 2+ )t' 2 (s 2A 



pc(h/4) = -f^ 2 - — ^ + - ^ 3 — — • (35) 

t 2 {s 2 +y \(s - s 2+ y {z -s 2+ yj \t' 2 {s 2+ y t' 2 ( S2+ y j\s-s 2+ zo-s 2+ ; 



Finally, we get 



with 



[ j " t 2 (s)-U(s)+A*oiAM {s y W 



A z«iam {s) = t 2 {s) (_ PC+ (1A 2 ) + PC{U/tl) + PC{l/t)) , (37) 
where, again, we evaluate l/t'(sA) in PC (1/t) perturbatively, obtaining for A(s) 

A z ° IAM (s) = A mIAM (s) - h^A mIAM (z a ) (38) 

t2{zny 

Formally, this looks the same as Eq. (|2"Tj) but now A mIAM (s) is of a more general form. Indeed, we obtain 
again A mIAM (s) plus a term depending on zq, which is also 0(/ -6 ). Nevertheless we will see that, as long 
as zq remains in the range of validity of the approximations, the latter is numerically small not only in the 
physical region but also in the subthreshold region, so that it can be neglected to obtain the naive simpler 
formula, thus justifying numerically the mlAM and the IAM. 

However, a derivation of the mlAM that differs from the standard one since subtractions are made at 
the Adler zeros is also possible for the unequal mass case, and we detail it next. 
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5.3 Dispersive derivation subtracting at the Adler zeros, 7 = 1/2, J = channel 

The derivation follows closely that in section 14.21 but now we will have extra terms (not negligible in the 
chiral expansion) coming from Discl/i2 hi the left cut when s A is expanded around s 2 . In addition, the 
expressions for the pole contributions are more complicated because i 2 '( s ) ^ 0. 
We again expand the left cut to NLO: 

LC (l/t) ~ 5 ' 52+ ' 54 / dz DiSC {1MZ \- WlM^l ( i + (39) 
2m J LC {z-s)(z-S2+) V Z ~ s 2+J 



-LC(n/t 



2\ 



S - S 2+ t 4 (S2-l 



t'(s 2 -)(s - S2-)(S2- - S2+) t' 2 (s2+)t' 2 {s2-){s2+ - S 2 -) 2 

where we obtain —LC(t^/t 2 ) plus two extra terms. The pole contributions now read 



I t"{s A ) 

t'(s A ){s-s A ) 2t'(s A y 



PC(Vt) = 7rrrxr „ — -^77^. (4°) 



pr ,, ,.2% _ *4pS2+) / ^( s 2+) t 4 (s 2 +)4 / ( g 2+) ^ 1 

t^Ly{l4 i 2 — -77 \2? Y2 + 777 \2 777 \3 y^ 1 

t2{S2+r{s- S 2 +) A \t2{S2+r * 2 ( s 2+)^ ) S - S 2 + 



+ 



t2(fl2+) t' A { S 2 + )t l 2 l (s2 + ) , 3t 4 (s2+)t 2 '(s2+) 2 t4(S2+)i 2 "(S2+) 



2^( S2+ )2 t 2 ( S2+ )3 4i' 2 ( S2+ )4 3t' 2 ( S2+ ) 3 



where all terms in PC(l/t) will be evaluated perturbatively excepting l/(s — s A ), that is the one that gives 
the needed pole at s A to the inverse amplitude. Then, adding the dispersion relations for 1/t and t^/t^, and 
using that 

1 = 1 , S-S 2+ t%{s 2 +) (42) 

t 2 (s) t' 2 ( S2+ )(s-S 2+ ) t' 2 ( S2 -)(s - S 2 -)(S2- - S 2+ ) 2t' 2 ( S2+ ) 2 [ ' 

we thus arrive to A mIAM in Eq. pTj) plus 0(/ -6 ) terms that can be safely neglected since they all correspond 
to the chiral expansion at very low energies. 

6 Numerical Results: comparison between different approaches 
and zo sensitivity. 

In this section we compare numerically the IAM with the modified methods we have derived. The precise 
values of the ChPT low energy constants are not relevant for what we want to show here. Just for illustration, 
we taken for SU(2) ChPT the typical values of 1% = 0.82, l\ — 6.2 from the second reference in [T], at a 
renormalization scale fi — 770 MeV. The values l\ = —3.7, l 2 = 5.0 have been obtained from an IAM fit 
to irir scattering data. In particular, we have used the same data sets used in the IAM fits of the second 
reference in [5] but we have updated the controversial / (600) channel by using the choice of data explained 
in [T^] consistent with forward dispersion relations and Roy equations as shown in [5U]. For wK scattering 
we have chosen the central values of the SU(3) ChPT constants of the IAM I set in the last reference of [5]. 

We show in FigfJ] several plots comparing the IAM and mlAM results for the modulus of different partial 
waves of irir and7r K elastic scattering in the scalar channel. We see that both methods are indistinguishable 
in the physical region, shown in the left column, and only differ in the Adler zero region which is shown in 
detail in the right column. Note that the IAM has a spurious pole and a double zero in that region, whereas 
the mlAM does not have such pole and its zero is simple. 

We have also calculated the a and k pole positions with the IAM and the mlAM, obtaining the same 
pole position within ~ 1 MeV, as shown in Table [1] 

In the derivation of the zoIAM, we have an arbitrary subtraction point zq ^ s 2 ,s A , but we only know how 
to calculate the amplitude at zq if we can use the chiral expansion l/t(zo) — 1/22 (zo) — t4,(zo)/t2(zo) 2 , which 
is only valid if zq lies on the low energy region. Also, due to t2 having a zero at s 2 , the above expansion is a 
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Figure 2: Comparison between the IAM and the mlAM for different isospin / partial waves of ww and irK 
in the scalar J — channel. The left column covers the region from the left cut up to 1 GeV. The only 
significant differences between both methods occurs in the region around the Adler zeros of each partial 
wave, which is shown in detail in the right column. 
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Method 


a pole 


k pole 


IAM 
mlAM 
z IAM, z = 8 t h 


443.71-? 217.58 
443.68 -i 217.56 
443.82 -i 216.99 


724.2 -i 216.2 

725.3 — * 216.3 
727.7-^210.0 



Table 1: a and k pole positions calculated with the IAM, the mlAM and the zqIAM with zq — s t h 



very bad approximation if zq is near s 2 . We can estimate how close zq could be from sa and s 2 by looking 
at the expansion 

1 ~ 1 U{Za) ~ 1 (\+ f 54 ,+■■■). (43) 



t(z ) t 2 (z ) t 2 (z ) 2 t 2 (z ) \~ ' {z ~s 2 ) 

where we have only made explicit the S4/ (zo — s 2 ) pole term. Hence, for our approximations to remain valid, 
we have to make sure that the ratio 54/(20 — s 2 ) is small enough. Thus, we show in Figl3]the contour plots 
in the energy and zq plane of the relative difference between the mlAM and the zoIAM 

\t^ AM {s)~t^ AM {s)\ 

\\t miAM {s)+t** iAM [s)y 1 ' 

as a function of the energy and zq. We show two contour lines corresponding to A = 10% and 5%. We see 
that, as long as the choice of z$ is sufficiently far from sa and s 2 (the white lines in the plots) the result 
of the zoIAM differs little from the mlAM. In particular, we have checked that in order to obtain a relative 
difference less than 10% and 5%, in the worst case, which is the J = 2, J = wave, we have to be sure 
that our choice of subtraction point zq makes the ratio of \s±/{zq — s 2 )\ < 3% and 2%, respectively. Let us 
recall that an error of 5% is a very precise result in this context, since it is much less than the uncertainties 
(mostly of systematic origin) of the existing data on meson-meson scattering. For instance, just the isospin 
violation effects, that are usually not taken into account by these experiments, can be estimated at the level 
of 2 or 3%, and these are added to further statistical and very big systematic uncertainties. 

In summary, the choice of subtraction constant has little relevance for the zoIAM, as long as it lies in a 
place where the NLO chiral expansion of the inverse amplitude is a good approximation to the inverse itself. 
In such case, the zoIAM results are very close to those of the mlAM, and therefore to those of the standard 
IAM in the physical region. This provides is a strong check of the stability and robustness of the standard 
NLO IAM results. 

We have also shown in the plots of Fig. 3 a line at z — s t h, which is a very natural choice of subtraction 
point, since the subtraction constants can then be written in terms of threshold parameters, which are well 
studied in the literature. We see that, with this choice of zq, in all cases we have a relative difference A 
which is less than 5%. Actually, it is even smaller, as seen in FigJH where we plot in detail the results for A 
subtracting at zq = s t h, and we see that for nn we have a relative difference less than 1.5% for all energies 
and less than 4% for the ttK case. 



7 Discussion and summary 

In this work we have shown that it is possible to modify slightly the one-channel Inverse Amplitude Method 
(IAM) so that, also in the scalar channel, it provides a reliable description of the unitarized partial wave 
amplitudes below threshold. In particular, we have shown that it is possible to obtain a modified IAM that 
has the Adler zeros located in the same place as the Effective Chiral Expansion up to the desired order, that 
these zeros are single, and that the spurious poles below threshold that occur in the standard IAM are no 
longer present. 

The IAM has been most frequently used at NLO, where the chiral expansion of a partial wave is written 
as t(s) — t 2 (s) + t4(s) + .... For such a case the simplest modification to the IAM that we have found can be 
written as follows: 

t miAM 1 \ MX) / 45 s 

[ ) t 2 { s )-t A { s ) + A™ iAM { s y [ ' 
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Figure 3: Contour plots of the relative differences between the mlAM and the zoIAM in the (s, zq) plane. 
We see that, as long as zq lies in the low energy region but sufficiently far from the the Adler zeros sa and 
s 2 (white lines) the differences become small for all partial waves. Note that for z = s t h (black line) the 
differences for all cases are less than 5%. 
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Figure 4: Relative differences A, with zo — sth, for nn (left), and 7rif scattering (right). Note that the 
differences are less than 1.5% and 4%, respectively. 



where 



\mIAM { „\ _ ^( S Y 



'00 



<2(*2) 



*4(S2) (S2-SA) (.,, s .,, , . t4(S2)t 2 '(«2) 

t 2 (s 2 ) - t 4 ( s 2) + 



(46) 



\S-S2? (S-S 2 )(S-S A ) \" " t' 2 ( S2 ) 

and in order to have the Adler zero exactly on its NLO position, we set sa — > s 2 + S4 where s 2 = 0{m 2 ) and 
s 4 ~ 0(m 4 /f 2 ) are the Adler zeros at LO and its NLO correction, respectively. That is, they are obtained 
from t 2 (s 2 ) = and i 2 (s 2 + S4) + t 4 (s 2 + S4) = 0. In general, S4 should be calculated numerically The above 
formula is valid both for the elastic scattering of equal and unequal meson masses. However, in most of the 
cases, like tttt scattering, the formula simplifies further since £ 2 (s 2 ) = and t 2 {s) = t 2 (s 2 )(s — s 2 ). Note also 
that the A mIAM (s) piece counts as next to next to leading order (NNLO) in the chiral expansion, and for 
that reason it was neglected in the standard IAM derivation, which is recovered by setting A mIAM (s) — ► 
and remains valid in the physical region, since there A mIAM (s) is indeed negligible. 

In Sect. 3 we have given a "naive" formal derivation of A mIAM (s) by adding, in a rather ad hoc way, 
the pieces needed to fix the Adler zero and spurious pole problems without spoiling unitarity and the chiral 
symmetry expansion. However, in Sects. 4 and 5 we have shown that the mlAM formulae above can be 
derived by using the analytic properties of amplitudes in the form of dispersion relations and imposing 
clastic unitarity on the right cut. The use of the Chiral Effective Expansion is well justified to calculate the 
subtraction constants and pole contributions to the dispersion relations, and is also used to approximate 
contributions from other cuts. In particular the integral over the left cut is calculated to NLO, which is a 
good approximation in the low energy region that dominates the integral. Therefore, there are no model 
dependencies but just approximations within the effective theory up to a given order. This allows for a 
straightforward and systematic extension of the elastic IAM and modified IAM to higher orders. 

A usual criticism to unitarization methods is their arbitrariness, but we have shown here that the IAM, 
modified or not, is not just unitarizing, but also imposing a stringent analytic structure on the amplitudes, 
something that leaves little room for such arbitrariness: the choice of the subtraction points. However, if the 
chiral effective expansion is to be used to calculate the amplitude at the subtraction points or pole positions, 
the subtraction points should lie, first of all, in the low energy region. But, in addition, since a dispersion 
relation is written for the inverse amplitude, that subtraction point should also lie far from the Adler zero. 
In this work, we have explicitly shown that as long as those two constraints are satisfied, the choice of 
subtraction point has a very small numerical effect on the resulting amplitude. Moreover, we have shown 
that the results of the standard IAM in the physical axis and resonance region in the complex plane remain 
unchanged when using the modified IAM and are extremely stable under different choices of subtraction 
points. 

In summary, we have presented a slightly modified Inverse Amplitude Method for the elastic case, that 
has the Adler zeros in the correct position and of the correct order and no spurious poles in that region. 
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We have shown that the results already obtained with the standard IAM are robust in the physical region, 
where it can still be used safely, but the new modifications allow for the study of the subthreshold region 
that is of interest in problems like thermal restoration of chiral symmetry [llj or the quark mass dependence 
[16j of resonances in meson-meson scattering. 
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A Chiral counting of pole contributions 

Next we will check by explicit calculation that the total pole contribution counts as 0(f~ 6 ) in the amplitude. 
We will do it for the more general case where t 2 (s) ^ 0. In particular, we only need to expand PC(\/t). 
First note that 

t'(s A ) = t' 2 {s 2 + Si + ...)+ t' 4 (s 2 + Si + ...) + 0(.r 6 ) = t' 2 (s 2 ) + 4(82)84. + t' A (s 2 ) + 0(f- 6 ). (47) 
With this 

1 _ 1 t'M , t'>(s 2 )t 4 (s 2 ) +oif _ 2y m 



t'(8 A ) t' 2 (8 2 ) t' 2 { S2 y t' 2 {8 2 f 

We also have to expand l/(s — s A ): 

1 - 1 ' 54 + 0(r 4 ) = — rn , 9 +0(f-% (49) 



s - sa s — s 2 (s - s 2 ) 2 ' s — s 2 t 2 (s 2 ){s - s 2 y 

and similarly for l/(zo — sa), where we have taken into account that t' 2 (s 2 )s4 = —t 4 (s 2 ) + 0(f~ 6 ). Hence, 
we find that 

Pcm -J-J-± l -) = -4 + -L-) - M (r-S - r-S) < 50 » 

t'(s A )\s-s A s-s 2 J t' 2 (s 2 )\s~s 2 20-S2/ t' 2 (s 2 y \(s - s 2 ) 2 {z -s 2 yj 

t>i{S2) < 2(s2)t4( S2 )^ ( 1 1 \ +0if -2 ) = p C{1/t2) _ pC{h/t 2 ) + 0{f -2 )i 



t^i 3 ?) 2 i 2( S 2) 3 /\ S ~ S 2 ZO ~ S 2 

so the total pole contribution PC{l/t)-PC{\/t 2 ) + PC{U/tl) = 0(/ -2 ), which yields a 0(/~ 6 ) contribution 
to t(s). 
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